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Abstract 

As the penetration level of transmission-scale time- 
intermittent renewable generation resources increases, 
control of flexible resources will become important to 
mitigating the fluctuations due to these new renewable 
resources. Flexible resources may include new or exist- 
ing synchronous generators as well as new energy stor- 
age devices. The addition of energy storage, if needed, 
should be done optimally to minimize the integration 
cost of renewable resources, however, optimal place- 
ment and sizing of energy storage is a difficult optimiza- 
tion problem. The fldelity of such results may be ques- 
tionable because optimal planning procedures typically 
do not consider the effect of the time dynamics of op- 
erations and controls. Here, we use an optimal energy 
storage control algorithm to develop a heuristic proce- 
dure for energy storage placement and sizing. We gener- 
ate many instances of intermittent generation time pro- 
flies and allow the control algorithm access to unlimited 
amounts of storage, both energy and power, at all nodes. 
Based on the activity of the storage at each node, we re- 
strict the number of storage node in a staged procedure 
seeking the minimum number of storage nodes and total 
network storage that can still mitigate the effects of re- 
newable fluctuations on network constraints. The qual- 
ity of the heuristic is explored by comparing our results 
to seemingly "intuitive" placements of storage. 



1 Introduction 

Electrical grid planning has traditionally taken two dif- 
ferent forms; operational planning and expansion or up- 
grade planning. The first is concerned with the rela- 
tively short time horizon of day-ahead unit commitment 
or hour-ahead or five-minute economic dispatch. The 
focus is on controlling assets that are already present on 



the system to serve loads at minimum cost while oper- 
ating the system securely. The second typically looks 
out many years or decades and is focused on optimal 
addition of new assets, again with a focus on mini- 
mizing the cost of electricity. When a system consists 
entirely of controllable generation and well-forecasted 
loads, the network power flows do not deviate signifi- 
cantly or rapidly from well-predicted patterns. In this 
case, expansion planning can be reasonably well sep- 
arated from operational planning. In the simplest ap- 
proach, expansion plans may be optimized against two 
extreme cases, e.g. the system's maximum and mini- 
mum load configurations. 

As the penetration of time-intermittent renewables in- 
creases, expansion and operational planning will neces- 
sarily become more coupled. For an electrical grid with 
large spatial extent, renewable generation fluctuations 
(here, we focus on wind generation) at well-separated 
sites will be uncorrected on short time scales[l , 2|, and 
the intermittency of this new non-controllable genera- 
tion will cause the patterns of power flow to change on 
much faster time scales than before, and in unpredictable 
ways. New equipment (e.g. combustion turbines or en- 
ergy storage) and control systems may have to be in- 
stalled to mitigate the network effects of renewable gen- 
eration fluctuations to maintain generation-load balance. 
It is at this point where operations planning must inter- 
face with expansion planning. The optimal placement 
and sizing of the new equipment, if required, depends on 
how the rest of the network and its controllable compo- 
nents respond to the fluctuations of the renewable gen- 
eration. Overall, we desire to install a minimum of new 
equipment by placing it at network nodes that afford us 
a high degree of controllability, i.e. nodes where con- 
trolled power injection and/or consumption have a sig- 
nificant impact on the network congestion introduced by 
the renewable fluctuations. From the outset, it is not 
clear which nodes provide the best controllability. Im- 



plicitly, placing a minimum of new equipment implies 
that it will experience a high duty, thus avoiding out- 
comes where equipment is only used for a small fraction 
of time. 

Before discussing our initial approach at integrating op- 
erational planning and expansion planning, we first sum- 
marize a few methods for mitigating the intermittency 
of renewable generation. When renewable penetration 
is relatively low and the additional net-load fluctuations 
are comparable to existing load fluctuations, a power 
system may continue to operate "as usual" with primary 
and secondary regulation reserves [3J being controlled 
via a combination of distributed local control, i.e. fre- 
quency droop, and centralized control, i.e. automatic 
generation control (AGC). In this case, planning for re- 
newables may simply entail increasing the level of re- 
serves to guard against the largest expected fluctuation 
in aggregate renewable output. 

As the penetration level grows, system operators may 
simply continue to increase the level of reserves, how- 
ever, if this fast-moving generation comprises natural 
gas combustion turbines, this simple planning will gen- 
erally result in increased renewable integration costs|4| 
which are usually spread over the rate base. Alterna- 
tively, operational planning can be improved by incor- 
porating long-range and short-range renewable genera- 
tion forecasting to better schedule the controllable gen- 
eration (energy and reserves) to meet net load and op- 
erate reliably in [5] El . Long-range forecasts, typically 
day-ahead or farther, are used in a unit commitment op- 
timization to ensure adequate generation and reserves 
will be online to meet the expected net load and its un- 
certainty. Short-range forecasts, typically hour ahead 
or shorter, are used in an economic dispatch optimiza- 
tion that sets actual generation and reserve levels of the 
committed generation. Some have investigated rolling 
unit commitmentllTlIll where updated wind forecasts are 
used modify the unit commitment more frequently. Sim- 
ulations have resulted in lower overall renewable inte- 
gration costs. 

Both unit commitment and economic dispatch seek min- 
imize the cost of electricity, however, they must also 
respect system constraints including generation limits, 
transmission line thermal limits, voltage limits, system 
stability constraints, and N-1 contingencies. Previous 
works iH |5] [71 |6l have generally looked at the effects of 
stochastic generation on the economics and adequacy of 
aggregate reserves while not considering such network 
constraints. These constraints may be respected for a 
dispatch based on a mean renewable forecast, however, 
if the number of renewable generation sites and their 
contribution to the overall generation is significant, veri- 
fying the system security of all probable renewable fluc- 



tuations (and the response of the rest of the system) via 
enumeration will be exponentially complex in the num- 
ber of sites resulting in a computationally intractable 
problem. 

The approaches summarized above do not consider net- 
work constraints or the behavior of the system on time 
scales shorter than the time between economic dis- 
patches (one hour in the case of [4]). In particular, they 
do not model how fast changes in renewable generation 
and the compensating response of regulation reserves in- 
teract with network constraints. In our initial study, we 
augment the approaches summarized above by focus- 
ing on the behavior of the electrical network at a finer 
time resolution and investigate how the control of energy 
storage affects its placement and sizing. We presume 
that the unit commitment problem has been solved, and 
at the start of a time period, we dispatch the control- 
lable generation by solving a DC optimal power flow 
(DCOPF) based on the site-specific mean forecast for 
wind generation. In the time before the next DCOPF 
is executed, the wind generation fluctuates, and we con- 
trol a combination of existing generators (using a simpli- 
fied description of frequency droop control) and energy 
storage to avoid violations of network constraints and 
generation limits. For each level of wind penetration, 
we generate many different realizations of wind fluc- 
tuations, and we gather statistics on the activity of the 
energy storage at each node. The statistics from simu- 
lated system operations are then coupled to the expan- 
sion planning process by developing a heuristic to guide 
the optimal placement and sizing of storage throughout 
the network-a result that cannot be achieved with the 
previous approaches described above. 

The rest of this manuscript is organized as follows. Sec- 
tion|2]lays out the system model including the develop- 
ment of our heuristic for placement and sizing of energy 
storage. Section [3] describes the simulations carried out 
on a slightly modified version of RTS-96l[8l and the ap- 
plication of our design heuristic. Section |4] gives some 
interpretations of the results, and Sectioii|5] wraps up 
with some conclusions and directions for future work. 



2 Mathematical Formulation 



We briefly describe the mathematical formulation of the 
optimal control problem and the algorithm used to solve 
it. We begin with a quick review of the DC power flow 
model and the proportional control scheme used to em- 
ulate frequency droop and AGCQ. 



2.1 Background and Notation 

Let Q = {'E,V) denote the undirected graph underlying 
the power grid, and n denote the number of nodes in the 
grid. Nodes are classified into three types: 

1 Loads t7/(assumed fixed over time scale of the 
study relative to the renewable generation). 

2 Traditional Generators Qg (whose output can be 
controlled). 

3 Renewable generators Qr (generators based on 
renewable sources whose output fluctuates over 
time). 

Given the vector p of power generated/consumed at ev- 
ery node in the grid, the power flow equations determine 
how much power f^ flows through each line in the net- 
work from node / to node /. For any 5 C 1^, we denote 
the power generation at nodes in 5 by p^ = {p, : / e 5} 
and use similar notation for other quantities(like en- 
ergy). The DC power flow equations are a linearization 
of the exact AC power flow equations, and have been 
widely used to compute optimal economic dispatch of 
generators in the grid. If Neb(y) denotes the set of 
neighbors of node j in the grid and x,j denote the induc- 
tance of the line from / to j, the DC power flow equations 
are given by: 

E9 , — 9, 



;GNeb(;) 



f,; = x,-,(e,-e,), 



where 0, is the phase of the complex voltage at node /. 
Written in matrix form, the first set of equations become 
p = — Lx9 where Lx is the Laplacian of the graph Q 
with edge weights l/x,,-. The Laplacian is not invert- 
ible, but if we are given a balanced power configura- 
tion, i.e., ^,p, — 0, we can always find a unique solution 
for the 9. This mapping can be constructed from the 
eigenvalue decomposition of Lx = ^/ w,e,ef by invert- 
ing the non-zero eigenvalues to get Mx = £,:„,. >o —^i^J 
and 9 — Mxp. Letting Mx' denote the /-th row of Mx, we 
can summarize the power flow equations as: 



f = Mp 

M'>=x,XMx'-Mx') 



(1) 



The aim of our energy storage control scheme is to dis- 
patch power from storage such that, in spite of fluctu- 
ations in the renewable generation, the following con- 



straints are maintained all the time: 

— f < Mp < f (Transmission Capacities) 

< Pt;' < ?§• (Generation Capacities) 
^p, =0 (Power Balance) (2) 

where the sum over / in the third constraint is a sum over 
three types of nodes. 

2.2 Power Generation and Proportional 
Control 

We will design a control scheme that operates in the du- 
ration T between two successive DCOPFs. At the begin- 
ning of each such time horizon, we assume that a subset 
of active generators Qg are committed (based on the pre- 
determined unit-commitment). We set their generation 
levels according to a DCOPF: 

min Y w;Pi 

Subject to Eqs.[T|andl2] 

where w,- are generation costs for each generator Let p" 
be the complete vector of initial generation at every node 
based on the DCOPF. The output p^^ of renewabnle gen- 
eration Qr is taken to be the mean of the renewable gen- 
eration p?, over the time between DCOPFs. We assume 
we know this value exactly, which implies we have per- 
fect predictions of the mean output of the renewables. 
In reality, the renewable generation P(^,. fluctuates dur- 
ing this time causing a mismatch between total genera- 
tion and total load. We assume that the generators Qg 
respond by changing their output by dividing the mis- 
match in some predetermined proportion a (usually cho- 
sen based on generation capacities)|9|. For notational 
convenience, we take a to be a vector of length «, with 
entries a, = for / ^ Qg. Let p'' = ^g^. — p^ denote the 
vector of renewable generation change. We assume for 
simplicity that loads are fixed, which is reasonable since 
in many cases fluctuations in renewables will be much 
larger than fluctuations in loads. The generators Qg will 
then respond as: 
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where non-zero components of ajy associated with gen- 
erators are chosen such that the net balance of power is 
maintained, l^p'' = "LieQgVi- 

2.3 Controlling Storage and Constraints 
on Storage 

We use s to denote the vector of energies stored at each 
node. However, the quantity we control directly is p*. 



the power drawn from each storage node at any given 
time. In our storage placement and sizing heuristic, we 
allow p* and s to be as high as needed to ensure that 
the constraints Q are satisfied. We also restrict the set 
of nodes that have storage to a set Qs- However, again 
for notational convenience, p*, s are still assumed to be 
vectors of length n and we set p^ = 0, S; = for j ^ Qs- 
In the paper, we use the optimal energy storage control 
as a heuristic to decide placement and sizing of storage. 
This decision only assumes that we know what the re- 
newable generation is going to be over the next win- 
dow of time of length T until the next DCOPF. Thus, 
we would like this decision to be robust to the behav- 
ior of wind beyond T and that our control over the next 
time window can work independently of what the con- 
trol was in the current time window. Hence we enforce 
the following invariance condition: No net energy is 
exchanged between the grid and the storage during the 
window of time T . Mathematically, this means that we 
require: 



s(o)-s(r) 



(3) 



We also assume that the storage is not accounted for in 
the DCOPF, i.ep^(0)=0. 

2.4 Overall Power generation 

If storage is being dispatched at a node, it adds to the 
pre-existing generation p"* such that total generation is 
p = p"* + p*, where p"' is the generation without in- 
cluding the local storage dispatch. For Qi or Qr, p'" is 
just the fixed load value or the time-variable renewable 
generation. The storage dispatch simply adds to these 
outputs. 

For generators, i.e. Qg, the situation is slightly more 
complicated. The dispatch of storage modifies the in- 
stantaneous load-generation imbalance so that the pro- 
portional control we proposed for Qg is slightly modi- 
fied. The controlled generation at a node responds to the 
sum of storage dispatch at all nodes, including its own, 
i.e. 

Thus, the overall power generation from a Qg node is 



p = p 



« -Ep^-EpH (4) 



2.5 



Optimal Control with perfect Fore- 
casts 





Figure 1 : Penalty Functions 



trol algorithm that assumes perfect forecasts of renew- 
able generation, i.e. the output of each renewable gen- 
eration node is known exactly over the time horizon T . 
The control algorithm uses this information to decide 
how much power to draw from or inject into each stor- 
age node in order to keep the grid within the constraints 
(|2]i. We note that, at the same time, the generation Qg 
is responding according to Q. The control algorithm 
computes this optimal dispatch by minimizing a cost 
function consisting of three terms. The cost function 
tries to enforce the constraints (|2| in a soft manner: We 
do this because in some situations, it might be infeasible 
to maintain all the constraints strictly. In those cases, 
we allow violations to occur but the controller suffers a 
cost that blows up cubically with the magnitude of the 
constraint violation. Thus, the controller will do its best 
to keep the system within its constraints, but when it is 
infeasible to do so, it will allow the violations to occur 
To make the notation compact, we define two auxiliary 
functions / and h, which are plotted in Figure [T| and de- 
fined below. 



f{x,a,b) = 







(k/* (x — a))^ 



if a <x < b 
ifx>b 
if X < a 



In this initial work on operations-based placement and 
sizing of energy storage, we implement an optimal con- 



h{x) ~ log (cosh (k/,jc)) 

where Kf — 50, K/, = .001. 

The three costs in the penalty function are: 

1 cli (p): penalizes overloading of transmission lines: 

cli(p)-I(,,,)eE/(M'^P,-f,„f,7) 

2 cg(p"*): penalizes violation of capacity constraints 
of controllable generators: 

cg(p"n=Le5,/(pr,0,P,) 

3 csp(p*): penalizes the absolute value of power 
drawn from storage (increases very slowly): 

csp(p^')=L-/'(P;) 

The fact that h increases very slowly ensures that 



cli(p) ,cg(p"'^) dominate this cost by far, so that 
the controller will never choose to allow violations 
to occur in order to save on power capacity. 

We define the total cost as cost(p,p"*, p*) — cli(p) + 
eg (p"') + csp (p*). The best control strategy minimizes 
this cost function integrated over time. In order to lever- 
age practical numerical techniques, we discretize the 
time axis uniformly with time step A to get Tf = \ time 
steps and assume that our controls are constant over each 
time period A. We seek 



P"5,(0),....P"5,(7/-1),=0 

subject to 

s(r) = X:p^-(OA 

x=0 



s(7»=s(0) 
pK0=O i^Q, 



p(o-p''+p^'(o+P'"(o+a(-Ep;'(o-Ep?(o) 



Levenberg-Marquardt correction IIT2I is required to en- 
sure convergence as the Hessian can become numeri- 
cally singular We use software from [ISJ to solve the 
problem . 

2.6 Optimal Sizing and Placement of Stor- 
age 

Renewable Generation Profiles 




p'"(f)=p(f)-p-'(r) 



Figure 2: Three typical renewable generation profiles 
used in our study. The dashed lines are the time averages 
of the profiles. 



Our problem falls within the well-studied class of de- 
terministic, discrete-time optimal control problems and 
can be solved using standard nonlinear programming 
techniques |10|. Note that p*^ and p'' are known from 
the DCOPF and the perfect renewable forecasts, respec- 
tively. Thus, wecanexpressp(f),p*(f),p'"'(f) directly in 
terms of p^ (f) using the second, third and fourth con- 
straints above to get p(p'^ )iP*(P(j ):P"(Pfi )■ Further 

we can write pV. \t) — -^ ^ — - — , with the conven- 
tion that Sg^{Tf) — S(^j(0) (in order to satisfy the con- 
straint ([3]l). This leaves us with the following uncon- 
strained optimization problem: 



mm 



Icl.fp('^-'' + '|-'<-'-"' 



eg P 



csp 



SgAt + '^)-Sg.Xt 



Sg,{t + l)~Sg^it) 



(5) 



It can be shown that this is a convex optimization prob- 
lem fill that can be solved efficiently using Newton's 
method (since the Hessian has a sparse block structure 
arising from the pairwise interaction). In practice, a 



We seek to develop heuristics to decide how to place 
storage and size its energy and power capacity. The 
high-level pseudocode given in algorithm [T] which uses 
the optimal control algorithm given above, is our initial 
attempt at this heuristic. Let Q^ denote the set of nodes 
with non-zero storage and perf(£7s) denote the inverse of 
the total storage power capcity used by the optimal con- 
trol with access to unlimited storage only at nodes in (j.,. 
At every iteration of the outer loop, the set of active stor- 
age nodes (^j is chosen. For each inner loop iteration, 
we create a trial where a random profile for renewable 
generation is created at each renewable node (details of 
how this is done are in Section [3] and typical examples 
of three profiles are shown in Fig. |2]). For each trial of 
the renewable profile, the control is allowed access to an 
unlimited amount of storage power and energy at all the 
nodes in ^s (and storage at nodes not in (jj), and we 
solve the control problem for the optimal storage power 
dispatch p*^ (f). For each trial k and each node in ^s, 
we find A'^, the maximum of the absolute value of the 
storage power at that node, which is a measure of the 
size of the storage power required at that node by the 
optimal storage dispatch for a given trial k. The values 
of A* are averaged over many trials. 
Next, the set of storage nodes Q^ is reduced by pick- 
ing only those nodes / e £7, such that A, > Ymax(A), i.e. 



Algorithm 1 Optimal Sizing and Placement Algorithm 

Choose thresholds e, e' 

g, ^ 1 : n 

repeat 

for A: = 1 ^ A^ do 

Generate Random Time Series Profiles for the 
Renewables 

Solve ^ with storage only at nodes £75 for the 
given profiles to get optimal dispatch p** (f) 
A*^ -(-max, |p"*(f)| 
end for 

y -;— maxJY : {peif ({/ e g^ : A,- > Ymax(A)}) < 

perf(^,)+e}} 

g,^{ie ^s:Ai>Ymax(A)} 
until y < e' 



only the nodes that require, on average, a storage power 
capability within a fraction y of the node with the maxi- 
mum storage power capability. In essence, we are using 
our optimal control scheme to identify the nodes that 
exercise the greatest control over cost(p,p"*,p'). The 
parameter y determines the severity of the "cut". As y is 
increased towards one, fewer and fewer nodes remain in 
g^- As our controllability is reduced, we might expect 
the level of constraint violations to increase. However, 
the violations seem to remain almost constant and we in- 
stead observe that at a critical value of y, the total amount 
of storage power on the network, as required by the opti- 
mal control, suddenly increases. We set y just below this 
threshold thereby selecting the reduced set of nodes that 
performs nearly as well as allowing the optimal control 
access to all of nodes. 

Using the reduced set of nodes from above as our new 
starting point, the algorithm is repeated (i.e. we regener- 
ate A and find the new maximal reduction in the number 
of nodes). This outer loop repeats until we are not able to 
reduce the number of nodes any farther without observ- 
ing a significant increase in the total amount of storage 
power. 



3 Simulations 

We tested our optimal control and heuristic for storage 
placement and sizing on a modified versions of RTS- 
96[81. The grid is shown in Fig. [3] Our modification in- 
cludes the addition of three renewable generation nodes 
shown in blue. The capacities of the new lines connect- 
ing the renewables to their immediate neighbors are set 
higher than the capacity of the added renewable gen- 
eration, otherwise, these Unes would be overloaded in 
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Figure 3: Our modified version of RTS-96. The added 
renewables are blue, loads are yellow and controllable 
generators are green. 



nearly every trial. 

3.1 Simulation Setup 

We performed three sets of simulations (corresponding 
to outer iterations of algorithm fTJ, and for each setting, 
we perform A^ = 2000 trials. We define the renewable 
penetration value of a trial to be 

The penetration is the fraction of the load that is served 
by the renewable generation. For each trial, a random 
value uniformly distributed between and 50% is se- 
lected for the penetration level, and the mean values of 
the renewables are scaled so as to achieve these values. 
A zero-mean fluctuation is generated around the base 
value for each renewable: 



p;-w-E 



1 Jq sm{&k{tA + z) +(^k) di 



where the (^k are randomly chosen phases and (Ok are 
chosen such that the above sine waves are harmonics of 
the basic wave with time period T. Since, we discretize 

^1 ^- • fn^ sin(cOi.(lA+-c)+(j)i.)£/x ^, 

the time axis, we use — — ^ ^ ^ — '—'^ — , the average 
power in the interval fA to (f + 1)A. The final profile 
is generated by scaling the disturbances so they are 
comparable to the mean generation (which is typical 
for turbulent fluctuations of wind) and adding it to 
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Figure 4: Activity Histograms: Red lines mark thresh- 
olds used for the reduction in the storage node set. 



the mean: pf'(f) = p°(f) * (1 +p;(f)),/ e Qr. Each 
instance of the renewable generation profiles we create 
is statistically independent and uncorrelated with every 
other instance. We do not intend for a sequence of these 
instances to approximate an actual time series of wind 
generation. 

We start with the unrestricted case (where unlimited 
storage is available at all nodes) and compute the av- 
erage activity of each node over all A^ = 2000 trials. The 
histogram of activities is shown in Fig. [4] The optimal 
y for the first outer iteration of algorithmfTlis 0.34 and 
the resulting reduced set contains 10 nodes. When we 
tried to remove more nodes from the set, the average 
level of total network power and energy storage capac- 
ity showed a sudden increase indicating that further re- 
ductions based on this first stage of planning actually 
decreased the system performance. 
We repeat this process, but this time starting with the 10- 
node set identified from the first iteration. Simulating 
another A^ = 2000 trials yields the second histogram in 
figure p] When we restrict our optimal control to just 
these 10 nodes, their rank order changes relative to the 
rank order in the unrestricted case. It is this property that 
allows our heuristic to make further reductions in the 
number of nodes without significant impacts to overall 
performance. 

The optimal y for this second iteration is 0.85 leaving 
only two nodes for the third interation. The algorithm 
terminates here as the controllable set cannot be shrunk 
any further without significance performance degrada- 
tion. The 3 cuts (all nodes, 10 nodes, 2 nodes) obtained 




Figure 5: Our modified RTS-96 grid showing the three 
sets of nodes identified by our heuristic. The minimal set 
(two nodes) is shown in red. The ten node set includes 
the two red nodes and the eight addtional green nodes. 
The maximal set includes all of the nodes. 



by the algorithm are depicted on the grid in Fig. [5] 

To study how the energy and power capacity require- 
ments of storage depends on the level of renewable pen- 
etration, we first must quantify these measures in a con- 
sistent fashion. For a given renewable profile p'^(f), the 
optimal storage power dispatch is p^*(f)^ and s*(f) is 
the resulting energy in storage. We define the energy 
in the renewable fluctuations to be s'^(f) — £!^^qP'^(x)A, 
i.e. s'^(f) is the energy stored in a (hypothetical) battery 
that is connected directly to a renewable node and elim- 
inates all fluctuations about the mean renewable genera- 
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Figure 6: The level of constraint violations vs the pene- 
tration of renewable generation. 
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Figure 7: The normalized energy capacity of storage in 
the entire network vs the penetration of renewable gen- 
eration. 
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Figure 8: The normalized power capacity of storage in 
the entire network vs the penetration of renewable gen- 
eration. 



tion. We will only look at changes in energy stored over 
time, so s*(0) and 8*^(0) can be arbitrarily chosen. With 
these brief prelimiaries, the parameters relevant to our 
study are: 

Normalized Power Capacity: This quantifies the to- 
tal power capacity of the storage relative to the sum 
of maximal power fluctuations over the renewables: 



IjGsimax, |p^*(f) 



LG^,(max, ?;-(?) -min,p[(f)) 

Normalized Energy Capacity: This quantifies the 
total energy capacity of the storage relative to the 
sum of maximal energy fluctuations over the re- 
newables: 

Iigsi(max,s*(f)-min,s*(f)) 

L-6§.(max,s;'(f) -min,s;(f)) 



Constraint Violations: This is defined as the sum of 
all constraint violations (as defined in (|2])) averaged 
over time. 

These quantities are all defined for a given renewable 
configuration. In order to investigate their typical be- 
havior as a function of renewable penetration, we create 
a 10 equally-spaced intervals between and .5. For each 
bin /, we compute statistics (mean and standard devia- 
tion) of all three quantities over trials with a penetration 
value in the bin /. The means an standard deviations for 
the three quantities are shown as functions of penetra- 
tion in Figs. 8|7|6 



4 Interpretations 

First, we note that Fig.|6]demonstrates the effectiveness 
of node-specific control for mitigating the impact of re- 
newable generation fluctuations on the network. The up- 
per curve in the inset in Fig|6] shows the level of con- 
straint violations for a system that maintains generation- 
load balance purely via proportional control, i.e. only 
the controlled generators responding to a generation- 
load imbalance in proportion to their capacity. Without 
the node-specific control in our model of energy stor- 
age, the response of the generators is inflexible, and sys- 
tem begins to experience constraint violations for small 
values of renewable penetration. With the node-specific 
control allowed by our storage model. Figure l6] shows 
that the control can eliminate all constraint violations 
up to about 30% renewable penetration for all the node 
configurations we explored. We do not attribute this be- 
havior to the storage itself, rather, this appears to be a 
property of the flexibility of node-specific control. Al- 
though we did not explore this possibility, we believe 
that if we would have incorporated the controlled gen- 
erators ^g into our control algorithm, we would have 
seen similar performance up to a renewable penetration 
of 30% for a system without any storage. It is interest- 
ing to note that more detailed studies have found similar 
thresholds for renewable penetration|4|. 
Figure l5] shows the location of the nodes with storage 
in the reduced-node sets identified by our heuristic in 
algorithm[T] The two-node set is perhaps the most illus- 
trative of how optimal placement and sizing of storage 
in a network does yield itself to simple rules of thumb 
or intuition. Intuition may have led to the conclusion 
that storage should be placed at the site of the fluctu- 
ating renewables because it can effectively mitigate the 
fluctuating power flows without these flows ever being 
injected into the network. It seems plausible that this 
should lead to a minimal amount of storage and quality 



control, however, the results in Fig. l5] and l6] show oth- 
erwise. The nodes in the two-node set do not include 
either of renewable generation sites. In fact, even when 
expanded to the ten-node set, only one of the three re- 
newable sites is included. These results demonstrate that 
the nodes that have a high degree of control over conges- 
tion caused by renewable fluctuations are not necessarily 
the renewable nodes themselves. 

To further demonstrate this point, we performed an ad- 
ditional simulation where the reduced-note set included 
just the three renewable sites. Figure |6] shows that the 
quality of control over the constraint violations is the 
same for this handpicked case as in the two-, ten-, and 
unrestricted-node sets. Figures l7] and IS] help to distin- 
guish between these node sets. For renewable penetra- 
tions below about 30%, the network total power capacity 
and energy capacity of the storage is nearly the same for 
each case. Beyond 30% where we begin to see a few 
constraint violations in Fig. l6j the two-node set shows 
significantly lower total network power and energy ca- 
pacity compared to the ten or unrestricted-node cases. 
The handpicked case is comparable to the ten-node case. 
A lower network total capacity implies that the nodes 
identified in the two-node case have the highest degree 
of controllability over the network congestion caused by 
the fluctuating renewables. 

A possible reason for the two-node set showing better 
controllability is evident from Fig. E\ These two nodes 
sit at the end or the middle of crucial transmission lines 
that link two major regions in RTS-96. In reference 
to Fig. |5] we name these regions "upper" and "lower". 
We conjecture that by controlling the power injections 
at these nodes, our control scheme can effectively con- 
trol to which region the renewable generation is directed. 
Therefore, our control can direct the renewable fluctua- 
tions to the region that, at that point in time, has the 
greatest ability to assist in mitigation. In addition to sim- 
ply absorbing the fluctuations locally, the storage may 
be used to redirect the renewable fluctuations through- 
out the network to leverage other resources to assist in 
the control. 



5 Conclusions and Future Work 

In this initial study, we have developed and demon- 
strated a method for coupling operations with system 
expansion planning for the optimal placement and siz- 
ing of storage in a grid with a significant penetration of 
time-intermittent renewable generation. Operations are 
incorporated into the planning via simulations of an op- 
timal control scheme that uses perfect renewable gen- 
eration forecasts to dispatch energy storage to eliminate 



violation of network and generation constraints. By sim- 
ulating the system on a relatively short time scale, we 
build up significant statistics on the storage activity at all 
network nodes. We have developed a heuristic that uses 
these statistics to reduce the number of storage-enabled 
network nodes while maintaining or improving system 
performance. 

Somewhat unexpectedly, our method does not preferen- 
tially place energy storage at the nodes with renewable 
generation. Instead, our method apparently favors nodes 
at critical junctions between major subcomponents of 
the network. We have proposed that these nodes pro- 
vide for enhanced controllability because, in addition to 
simply buffering the fluctuations of the renewables, con- 
trolled power injections at these nodes can modify over- 
all network flows and direct fluctuating power flows to 
regions that are better positioned to mitigate them. 
There is much follow on work needed to expand the con- 
cept presented in this manuscript and to verify some of 
its conjectures. A few examples include: 

• Additional networks should be considered, includ- 
ing different configurations of renewable genera- 
tion on the present network, to investigate whether 
our node-reduction heuristic is robust. 

• The explanation of the placement of storage pro- 
posed in this work should be verified by correlat- 
ing the response of the storage with nearby power 
flows. 

• The operational simulation should be made more 
realistic by incorporating ramping constraints on 
the controllable generation. 

• The method should be extended to consider alterna- 
tives to proportional control for changing the out- 
puts of the controllable generators. 
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